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ABSTRACT 



This is a sequel to the previous paper, in which we investigated the structure 
and stabihty of the spherically symmetric accretion flows through the standing 
shock wave onto the proto-neutron star in the post-bounce phase of the collapse- 
driven supernova. Following the prescription in the previous paper, we assume 
that the accretion flow is in a steady state controlled by the neutrino luminosity 
and mass accretion rate that are kept constant. We obtain steady solutions for 
a wide range of neutrino luminosity and mass accretion rate. In so doing, as 
an extension to the previous models, we employ a realistic EOS and neutrino- 
heating rates. More importantly, we take into account the effect of convection 
phenomenologically. For each mass accretion rate, we find the critical neutrino 
luminosity, above which there exists no steady solution. These critical points 
are supposed to mark the onset of the shock revival. As the neutrino luminosity 
increases for a given mass accretion rate, there appears a convectively unstable 
region at some point before the critical value is reached. We introduce a phe- 
nomenological energy flux by convection so that the negative entropy gradient 
should be canceled out. We find that the convection lowers the critical neutrino 
luminosity substantially, which is in accord with the results of multi-dimensional 
numerical simulations done over the years. We also consider the effect of the 
self-gravity, which was neglected in the previous paper. It is found that the 
self-gravity is important only when the neutrino luminosity is high. The criti- 
cal luminosity, however, is little affected if the energy transport by convection is 
taken into account. 
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Subject headings: stars: — shock waves — hydrodynamics — supernovae: general 

1. Introduction 

Core-collapse supernovae are triggered by the gravitational collapse of massive stars. 
The mechanism of explosion is still not well understood. Many researchers think that the 
so-called neutrino heating mechanism is the most promising at present. In this scenario, 
after the shock is launched by the core bounce, it stalls inside the core due to the energy 
loss by photo-dissociations of nuclei as well as neutrino-cooling, and then it is revived by 
the irradiation of neutrinos diffusing out of the proto-neutron star (Wilson 1985). Although 
many researchers have studied this scenario, numerical simulations have shown that the 
stalled shocks do not revive as long as the collapse is spherically symmetric (Licbcndorfcr et 
al. 2001; Buras et al. 2003; Thompson et al. 2003; Liebendorfer et al. 2005; Sumiyoshi et al. 
2005). 

After the shock is stagnated due to the energy losses, the accretion through the standing 
shock wave onto the proto-neutron star is approximately steady Burrows & Goshy (1993) 
and Yamasaki & Yamada (2005) attempted to approximate this phase of evolution by so- 
lutions of the time-independent Euler equations, assuming constant mass-accretion rates 
and neutrino luminosities. Adopting such a method, they demonstrated clearly what is the 
cause of the failure of explosions found in the numerical simulations. Varying these con- 
trolling parameters, the mass-accretion rate and neutrino luminosity, they found that for a 
given mass-accretion rate there is a critical neutrino luminosity, above which there exists 
no steady solution. Based on this fact, they argued that the revival of stalled shock occurs 
when the neutrino luminosity exceeds this critical value. 

Yamasaki & Yamada (2005) found further that there are in general two types of solutions 
for a given mass accretion rate and a neutrino luminosity. The shock radii in these two 
solutions differ. They referred to the solution with a smaller shock radius as the inner solution 
and to the other solution as the outer solution. As the neutrino luminosity is increased with 
the mass accretion rate fixed, the shock radius for the inner solution becomes larger, whereas 
that for the outer solution becomes smaller. Two solutions coincide with each other when 
the luminosity reaches the critical value. For the luminosity over the critical value, there is 
no solution. They also studied the linear stability of the system against radial perturbations 
and found that the inner solution is stable while the outer solution is unstable and they are 
merged at the critical luminosity and become marginally stable. In that paper, the equation 
of state and the reaction rates for neutrino processes were simplified, and the self-gravity of 
the accreting matter is neglected. We improve these aspects in this paper and discuss their 
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effects on the results. 

The other concern in this paper is the convection in the accreting flow. The convection in 
the supernova core and its imphcation for the explosion has been studied by many researchers 
using multi-dimensional numerical simulations over the years (Herant et al. 1994; Burrows 
et al. 1995; Janka & Miiller 1996; Mezzacappa et al. 1998). They have demonstrated that 
the convection tends to assist the revival of shock owing to the non-radial motions and 
the enhanced heating of the accreting matter although it seems that the convection alone 
can not give a successful explosion. Our main concern here is whether we can understand 
in the frame work by Yamasaki & Yamada (2005) the results obtained by the large scale 
simulations . In fact, we have recognized that some solutions in the paper (Yamasaki & 
Yamada 2005) have a region that has a negative entropy-gradient and is unstable against 
convective motions. In this paper, we discuss the effect of convection by incorporating a 
phenomenological energy flux in our time-independent, spherically symmetric models. The 
convective energy flux is introduced in such a way that the negative entropy-gradient should 
be canceled out. This assumption will correspond to the maximal mixing by convection and, 
hence, we can hopefully estimate the upper limit for the effect of convection on the critical 
neutrino luminosity. 



2. Models 

In this section, we describe our models, giving the basic equations and some assumptions. 
We are interested in the post bounce phase, in which the shock is stalled inside the core and 
becomes a standing accretion shock. Matter is accreted onto the proto-neutron star through 
the shock wave, being irradiated by neutrinos diffusing out of the proto-neutron star. 

Assuming that this accretion flow is steady and spherically symmetric, we solve the 
following time-independent ordinary differential equations. 

47rrV«r = M, (1) 

dur I dp GM ^ 
dr p dr 

dS . . . 

'^rpT-j^ = q, (3) 

dYe , 

UrU—— — A, (4j 

dr 

— = 47rrV, (5) 
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where r, Ur, p, p, -S", Ye and n denote the radius, radial velocity, density, pressure, entropy 
per unit mass, electron fraction and baryon number density, respectively; M, G and M are 
the mass accretion rate, gravitational constant and mass enclosed in the sphere of radius r, 
respectively; q and A are the heating and reaction rates due to the reactions with neutrinos, 
respectively. The rotation and magnetic field are neglected. The Newtonian formulation is 
justified because the region of our interest is outside the proto-neutron star and the general 
relativistic correction is at most 10%. 

Wc do not consider the flow outside the shock wave, simply assuming a free fall, and 
restrict our calculations to the region inside the shock. We impose the Rankine-Hugoniot 
relations at the shock front, r = Tgi^: 

pul+p^ puul+p^, (6) 

-u^ + e+ - = -u^ + eu+ —, 7 
2 P 2 pu 

where e denotes the internal energy including the nuclear binding energy. The suffix u 
stands for the upstream quantities given as follows. Assuming that the upstream matter is 
composed of cold irons and flowing into the shock at a free-fall velocity, we take 

= 0, (8) 

and 



u^ = ^——. (10) 

The inner boundary of the calculation region is set at the neutrino sphere, from which, wc 
assume, thermal neutrinos of all flavors with a temperature of = 4.5MeV are radiated 
outwards. The radius of the neutrino-sphere, r,^, satisfies the following relation, 

L.^ ^ ^aT,' ■ ATrrl (11) 
where a is the Stefan-Boltzmann constant. At the inner boundary, we also assume 

M = I.SMq, (12) 

that is, the mass of the proto-neutron star. 

As for the neutrino reactions, we take into account only the emission and absorption by 
free nucleons. 

p + e~ ^ Ue + n, (13) 
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n + e+^i/e+P- (14) 
The heating and reaction rates can be decomposed accordingly as 

Q^Qep- <iun + Qe+n " Qpp: (15) 
A = —Aep + A^n + Ac+n ~ App. (16) 

These rates are calculated based on the formulae given by Bruenn (1985). In the previous 
paper (Yamasaki & Yamada 2005), a simplified equation of state was employed and the 
photo-dissociations of nuclei are neglected for simplicity. In this paper, we employ a realistic 
equation of state by Shen et al. (1998) and take into account the nuclear composition and 
dissociation consistently. 

The optical depth of the inner boundary (or the neutrino sphere) from the shock front 
is set to be 2/3 for the electron- type neutrino. 

Tu, = J a^Jr = -, (17) 

where a^,^ is the energy-averaged absorption coefficient for the reactions of the electron-type 
neutrino given in equation (13). We neglect the contribution from the layer outside the 
shock wave to the optical depth because the absorption coefficient is negligibly small there. 
Since the region of our interest is outside the neutrino sphere, we do not solve the neutrino 
transfer equations, assuming that the luminosity and energy of neutrinos are independent of 
radius. 

Adopting the neutrino luminosity, L^^, and the mass accretion rate, M, as model pa- 
rameters, we solve equations (l)-(5) together with the boundary conditions for a wide range 
of these parameters. The phenomenological treatment of convection will be described in 
section 3.2. 



3. Results 

3.1. Effects of Improved Physics 

First, we show the results of the calculations with the self-gravity neglected. We fix the 
mass to be M = 1.3M0 and ignore equation (5). The radial infall velocities of the spherically 
symmetric steady flows are given in flgure 1. Although we have obtained solutions for various 
values of mass accretion rate, we show only the case with M — 1.0MqS~^. For other values 
of M, the results are similar qualitatively. As in the previous paper (Yamasaki & Yamada 
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2005), we find that for a given mass accretion rate and neutrino luminosity, there are two 
types of solutions, each having different radii (that is, the inner solution and outer solution) 
as long as the neutrino luminosity is below a certain critical value. As shown in figure 1, 
the shock radius for the inner solutions becomes larger, whereas that for the outer solutions 
becomes smaller, as the neutrino luminosity is raised with the mass accretion rate fixed. 
These two solutions coincide with each other when the luminosity reaches the critical value, 
and there is no solution for the luminosity exceeding the critical value. 

This behavior of the solutions is understood as follows. Equations (l)-(4) can be re- 
assembled for the radial velocity as, 

i2_f9p\ ]d\n\ur\_2/'dp\ _GM_l{rdp\ q / 9p \ a1 
\dpJs,Ye) ^\9p)s,Ye P \\9S J p^Ye '^rpT \dYeJpsUrnf' 

(18) 

As known for the Bondi fiow (Shapiro and Teukolsky 1983), the subsonically accreting matter 
is first accelerated and then decelerated, giving a peak in the r — u,,. diagram (see the region 
with log^^r ~ 3.5 in figure 1). [Note that the post shock matter is not hydrostatic and the 
advection cannot be ignored particularly when the shock radius is larger than about lO^cm 
for the mass accretion rate of 1 Mq/s.] This is because the first term on the right-hand side 
of the equation (18) dominates the other terms in the outer region, whereas the second term 
becomes dominant in the inner region. The third and forth terms on the right-hand side are 
neghgible at large radius. Since the free-fall velocity decreases monotonically with r, there 
are in general two solutions satisfying the Rankine-Hugoniot relations. When the neutrino 
luminosity increases, the accretion velocity is lowered, and the solution ceases to exist at 
some point, which gives the critical neutrino luminosity. In the present case, however, the 
behavior is more complicated. The accretion velocity has a second peak (see the region with 
logio ~ 2 in figure 1) when the luminosity becomes large and the heating zone appears. The 
emergence of the second peak is attributed to the neutrino heating and cooling represented 
by the third and forth terms on the right-hand side of equation (18). Since {dp/dS)p^Ye 
is positive and the last term is always smaller than the third term on the right-hand side 
of equation (18), the non-adiabatic effects represented by the last two terms work so as to 
reduce the deceleration of fiows in the heating region. In addition, as the infiowing matter 
is heated, the value of {dp/dp)s,Ye becomes larger, and if the heating is great enough, the 
sum of the right-hand side of equation (18) becomes positive at some point and the fiow 
switches from deceleration to acceleration. After matter enters the cooling region, the value 
of {dp/dp)s,Y^ decreases again and the second term on the right-hand side of equation (18) 
dominates over the other terms at some radius and, as a result, the flow is decelerated again 
in the innermost region (the upper left panel of figure 2) . The improvement of the EOS and 
the proper treatment of the photo-dissociations of nuclei are responsible for the difference 
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from the results in the previous paper. The temperature is substantially smaller in the 
present models. Apart from this non-monotonic behavior, however, the characteristics of 
the solutions found in the previous paper still hold. We think that the outer solution will 
not be realized in the actual evolution of the supernova core because the shock radius for the 
outer solution is extremely large for rather small luminosities that are commonly found by 
detailed numerical simulations in the phase of the shock stagnation. Hence we will discuss 
only the inner solutions from now on. 

The results with self-gravity arc shown in figure 2. As long as the neutrino luminosity 
is less than about 12 x 10'''^crgs s^'^, the effect of self-gravity is negligible, and the solution 
curves are almost the same as those neglecting the self-gravity, while the effect becomes 
significant quickly when the luminosity exceeds this value. As can be seen in figure 3, 
where we show the shock radius and the mass enclosed in the shock surface, they increase 
rapidly as the luminosity approaches the critical value. The effective increase of gravitational 
attraction owing to the incorporation of self-gravity in the second term on the right-hand side 
of equation (18) tends to decrease the velocity of the subsonic flow. This may appear to be 
counter-intuitive, but can be understood as follows. When the accretion flow is supersonic, 
the pressure is less important and the increase of mass simply results in the increase of 
flow speed. For the subsonic flow, on the other hand, the pressure plays an important role 
and the increase of gravity results in the increase of pressure as well as the density, which 
in turn leads to the decrease of the flow velocity because of the mass conservation. The 
effect of self-gravity is negligible for the inner solutions in most cases, while it seriously 
affects the outer solutions, in which the flow is initially accelerated. Since the density is 
increased, the shock radius becomes smaller so that the optical depth for neutrinos be 2/3 
at the inner boundary (17). This in turn leads to the smaller critical luminosity, at which 
the outer solution coincides with the inner solution and the steady solution disappears. It 
is interesting to see that the shock radius for the critical luminosity is almost unchanged as 
the mass accretion rate is varied. It becomes slightly larger when the mass accretion rate 
gets smaller than 1.0 Mq. 

In flgure 4, we show the critical luminosities for various mass accretion rates, both with 
and without self-gravity taken into account. As shown also in flgure 2, there emerges a region 
with a negative entropy-gradient when the luminosity exceeds a certain value for each mass 
accretion rate. Such a region is convectively unstable by the classical criterion (see section 4 
for possible corrections to the criterion) , We will discuss the effect of convection in detail in 
section 3.2. Since the self-gravity becomes appreciable quickly once the luminosity exceeds 
a certain value that is slightly smaller than the critical luminosity, the critical luminosity 
is signiflcantly affected by the self-gravity as long as the convection is ignored. Comparing 
flgure 4 with the previous results (flgure 1 of Burrows & Goshy (1993) or flgure 8 of Yamasaki 
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& Yamada (2005)), we can see that the critical luminosity is substantially increased in the 
present models. In section 3.2 we will discuss how this will be changed again after the 
convection is taken into account. 



3.2. Effect of Convection 

In subsection 3.1, we have seen that the accretion flow has in general a region with a 
negative entropy-gradient if the neutrino luminosity is larger than a certain value. Then, we 
expect and multi-dimensional numerical simulations have demonstrated that the convcctive 
motions occur. The convection induces not only the out-going energy flow effectively but 
also the mixing of matter, and, as a result, affects the structure of the accretion flow and the 
location of the shock. In this section, we discuss the effect of the convection on the steady 
flow as well as on the critical neutrino luminosity in the present frame work. 

Recently, Foglizzo et al. (2006) pointed out that the negative entropy-gradient is not a 
sufficient condition for the convection in the accreting matter and that the convection will be 
suppressed by the advection. In section 3.1, we used a classical criterion to judge if there is a 
convection zone or not, which would have been an overestimation for the convectivc region. 
In this section, we still stick to the classical criterion, since the approach proposed below is 
highly phenomenological and the results we obtain should be regarded as " quahtative" , but 
not "quantitative". The possible suppression of convection due to advection in our models 
will be assessed in section 4. 

In order to incorporate the effect of convection in our spherically symmetric models, we 
introduce the effective energy- and proton-number fluxes as follows. In the convection zone 
{t^ > '"gain), the equations for energy (3) and electron fraction (4) are replaced by 

m ■ i. Ul 1 exonv /ir\\ 

UrpT— = g - — , (19 

dr dr 

dY 1 dr'^F 

"'-'e \ ± uji -I n.conv /oaN 

UrU— = A - — , (20) 

dr r^ dr 

where Fg^conv and F„ conv are the energy- and the proton-number fluxes associated with the 
non-spherical convective motions, respectively. 

It is admittedly difficult to give these fluxes from the first principle. Multi-dimensional 
numerical simulations done so far suggest that the angle-averaged treatment like the one 
employed in this paper may not be a very good approximation for the convection in the 
heating region. With these caveats in mind, however, we proceed as follows. Recalling that 
the convection is driven by the negative entropy- gradient and tends to reduce it, we assume 
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in the convective zone that the right-hand sides of equations (19) and (20) vanish, that is, 

. 1 dr^F,, 



J.2 ^j. 

and 



= 0, (21) 



1 dr'^F 

^_ 1 ar J-n,conv .22) 

This is in a sense the hmit of maximally efficient mixing, which realizes an isentropic distri- 
bution in the convective zone. We regard these equations as the definitions of the fluxes. 

Since the convective region is extended up to the shock front, these fluxes should be 
also taken into account in the jump conditions for the shock wave. We modify equations (7) 
and (9) in the original condition to 

1 2 I , P , -^e,conv ■'■ 2 , i ^'u ^r,o\ 

-u^ + e + - + ^ = -M^ + euH , (23) 

Ye + = i;,u. (24) 

nur 

We further impose the following conditions, 

-^e,conv -^n,conv 0- (25) 

at the inner boundary of the convective zone, that is, the gain radius r = rgain, where 

q = 0, (26) 

is satisfied. 

The solutions of the above equations are shown in figures 5 to 7, where the self-gravity 
is taken into account. In figTirc 5, wc can sec that the complicated post-shock velocity distri- 
butions with a doublc-pcak have been replaced by the monotonic ones after the convection 
is taken into account. This is understood as follows. In our treatment of convection, the 
last two terms on the right-hand side of equation (18) is just compensated by the convective 
flux in the heating region. As a result, the situation in the heating region becomes very 
much similar to the adiabatic flow, and the flow is monotonic again. More importantly, the 
critical luminosity is lowered by a factor of ~ 2 by the convection as demonstrated in flgure 
6. This agrees qualitatively with the results of multi-dimensional numerical simulations. In 
terms of our model, the reason is the following. The convection tends to raise the entropy 
just behind the shock wave. This leads to the decrease in the accretion velocity, which then 
makes it easier to reach the critical point. 
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It is also noted that the shock radius for the critical luminosity becomes smaller (see 
figure 7), compared with the models neglecting the convection. This implies that the results 
are hardly affected by self-gravity if the convection is taken into account, since the self- 
gravity becomes important only when the critical point is approached in the models without 
convection. 



4. Summciry and Discussion 

In this paper, we have obtained spherically symmetric steady accretion flows through 
the standing shock wave onto the proto-neutron star. We have improved the treatment of 
microphysics such as the equation of state and neutrino reaction rates. We have found that 
our previous findings arc essentially unchanged after these modifications. For a given mass 
accretion rate, there are in general two solutions, the inner and outer solutions. There is a 
critical neutrino luminosity for a given mass accretion rate, above which no steady solution 
exists. The outer solution is unstable against radial perturbations and the two branches 
of solution coincide with each other for the critical luminosity. There is a quantitative 
difference, though. The critical luminosity, in particular, becomes much higher when the 
realistic microphysics are implemented. 

We have also investigated the effect of convection, employing a phenomenological de- 
scription of convection. We have introduced the convective fluxes of energy and proton 
number, which then have been determined in such a way that the negative entropy gradient 
is canceled out by them. This is roughly corresponding to the maximally efficient mixing 
by convection. We have found that the critical luminosity is lowered by a factor of ~ 2 
by the convection, consistent with results from the large-scale numerical simulations done 
so far. The latest results of detailed numerical simulations suggest that the convection is 
not sufficient for the shock revival (Buras et al. 2006). Comparing the values of the criti- 
cal luminosity in our models with the actual luminosities obtained by detailed simulations 
(Sumiyoshi et al. 2005), we have found that the former is larger. Thus, our results agree 
with the numerical simulations also in this respect. It is noted that our models are supposed 
to correspond to the maximal convection. Furthermore, we have found that the profiles of 
various quantities, e.g. the density profile with as well as the shock radius obtained 

here are in good agreement with the results of detailed simulations for relatively low lumi- 
nosities of < 6 X lO^^ergs s^^. For higher luminosities, however, our results differ from those 
of the simulations, suggesting that the latter are far from the critical lines for shock revival. 

In this paper, we have employed the classical criterion for the convective stability, which 
is, rigorously speaking, applicable only to the conflguration without advection. Recently, 



- 11 - 



Foglizzo et al. (2006) gave the modified criterion that takes into account the advection. 
They introduced the parameter x defined as 



r 


N 


' ' gam 


Ur 



X = 



where N is the Brunt- Vaisala frequency given as 



dr 



dr, 



(27) 



dp 



dY, 



9 



P \9SJ p,Y. P p,s dr J Ti 

In the above expression, Fi and g are given by 

dlnp^ 



(28) 



and 



9 



dlnp 



GM 



S,Ye 



(29) 



(30) 



They claimed that the convection occurs in the presence of advection if the entropy-gradient 
is negative and x is larger than 3. In order to see the consequence of this modification to 
our models, we have calculated the values of x ioi the solutions without convection taken 
into account. The values of x for our models with self-gravity are shown in figure 8. We 
can see that x is always smaller than 3 for the mass accretion rates larger than I.OMqS"^ 
and that even for the mass accretion rate of O.IMqS"^, x is smaller than 3 except for a 
narrow range of neutrino luminosity. Thus, taken at face values, the criterion by Foglizzo et 
al. (2006) predicts that the convection does not occur for most of our models. This seems 
to be at odds with the results of numerical simulations. It is, we think, worth notice that 
the steady solutions we have obtained here are qualitatively different from the unperturbed 
states assumed in Foglizzo et al. (2006). Foglizzo et al. (2006) assumed simplified formulae 
for the neutrino heating and equation of state while we employed more realistic ones and, 
as a result, the post-shock fiows are much complicated. Hence we suspect that the critical 
value of X may be lowered in this case. Indeed, the detailed linear analysis using our steady 
solutions suggests that the critical value of x may be smaller than 3 (Yamasaki & Yamada 
2006). However, the present models employing the classical criterion will give the upper 
limit for the effect of convection. 

There are a couple of things that have not been considered in this paper but will be 
important. For example, we have introduced the phenomenological energy- and proton- 
number fiuxes in the convection region, but have neglected the kinetic energy of turbulence. 
This may not be ignored for the large velocity fiuctuations as observed in the numerical 
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simulations. Taken into account, it might increase the critical luminosity because the energy 
of down-flows will be sucked. It is also pointed out that the non-spherical deformation of 
the shock front by convection may be an important ingredient that cannot be considered in 
the present frame work. Thus far we have paid attention only to the condition for the shock 
revival. This is, however, a minimum requirement. The explosion energy and the mass 
of neutron star left behind are equally important issues. We are currently attempting to 
estimate these values based on the models presented here. Our model is also being extended 
to accommodate rotation, magnetic field, anisotropic neutrino irradiation as well as being 
utilized for the linear analysis of the non-local instability (Foghzzo 2002; Blondin et al. 2003; 
Ohnishi et al. 2006), whose results will be pubhshed in the near future. 

This work is partially supported by the Grant-in- Aid for the 21st century COE program 
" Holistic Research and Education Center for Physics of Self-organizing Systems" of Waseda 
University and for Scientific Research (14740166, 14079202) of the Ministry of Education, 
Science, Sports and Culture of Japan. 
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Fig. 1. — Left. Accretion velocity for the inner solutions with M = 1.0 Mq s ^, Lj,^ = 
(2,4,6,8, 10, 12, 14, 16, 18) x 10^^ ^^-gg g-i ^nd 18.7125 x 10^^ g^-gg g-i (critical value) from 
left to right. Right. Accretion velocity for the outer solutions with M — 1.0 Mq s~^, 
L^^ = (17.0, 17.5, 18.0, 18.5) x 10^^ ergs s"^ and 18.7125 x 10^^ ergs s"^ from top to bottom. 
Dashed lines denote shock jumps. Dotted lines show the downstream velocity satisfying the 
Rankine-Hugoniot relations at each radius. The self-gravity is not taken into account. 
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Fig. 2. — Accretion velocity (upper-left panel), density (lower-left panel), entropy 
(upper-right panel), electron fraction (lower-right panel) for M — 1.0 Mq s~^, L,y^ = 
(2, 4, 6, 8, 10, 12) X 10^2 gj-gs g-i ^nd 12.4606 x 10^^ gj-gg g-i (critical value) from left to 
right. 
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Fig. 3.— Left. Radius of the shock surface for M = 0.1, 1.0, 2.0, 3.0, 4.0 s"^ from left to 
right. Right. Mass enclosed in the shock surface for M = 0.1, 1.0,2.0,3.0,4.0 Mq from 
left to right. 
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Fig. 4. — Left. Critical luminosity (solid curve) and boundary for convective/non-convective 
solutions (dashed curve). Self-gravity is taken into account. Above the dashed curve, there 
exists a convectively unstable region in the solution. Right. The same as in a. but without 
self-gravity. 
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Fig. 5. — As in Fig. 2, but with convection taken into account and with M = 1.0 Mq s~^, 
Li,^ = (4,5,6) X 10^^ ergs s^^ and 6.7860 x 10^^ ergs s^^ (critical value) from left to right. 
Dashed lines denote the shock jumps. 
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Fig. 6. — Critical luminosity (solid curve) for the models with convection taken into account. 
Above the dashed curve, there is a convection zone in the solution. 
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Fig. 7. — Radius of the shock surface (sohd curves) for the models with convection taken 
into account. The mass accretion rates are M — 0.1,1.0,2.0,3.0,4.0 Mq from left to 
right. Dashed curves denote the gain radius. 



-20- 



X 




100 200 300 

L^^ [10^^ ergs/s] 



Fig. 8.— The parameter x for M = 0.1,1.0,2.0,3.0,4.0 M© s"^ from left to right. The 
self-gravity is taken into account for these models. 



